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Abstract-  Multipath  interference  is  a  major  source  of  noise 
for  synthetic  aperture  sonar  systems  operating  in  shallow 
water.  Motivated  by  this  problem,  we  present  an  iterative 
algorithm  for  blindly  separating  the  signal  from  the  multipath 
interference  that  uses  differences  in  the  temporal  coherence 
properties  of  the  signal  and  multipaths  caused  by  sea  surface 
roughness  to  estimate  the  optimum  filter  weights.  The  filter 
weights  are  estimated  by  minimizing  the  circular  variance  of 
the  phase  differences  between  overlapping  vertical  phase 
centers.  An  important  advantage  of  this  approach  is  that 
signal-free  training  data  and  accurate  array  calibration 
information  are  not  needed.  Experimentally  we  show  that  the 
blind  separation  performance  is  competitive  with  high 
resolution  angle  of  arrival  estimation  and  compares  favorably 
to  the  Cramer-Rao  lower  bound  predicted  error. 

I.  INTRODUCTION 

Synthetic  aperture  sonar  (SAS)  is  commonly  used  in 
applications  requiring  high  resolution  seafloor  images  such 
as  bottom  mapping,  mine  hunting,  cable  and  pipeline 
surveying,  and  underwater  archaeology.  Fine  cross-range 
resolution  is  obtained  using  the  motion  of  the  sonar  platform 
to  synthesize  a  much  larger  virtual  acoustic  aperture  or 
baseline  than  the  sonar’s  physical  aperture  by  coherently 
combining  echoes  from  multiple  pings  [1].  In  addition  to 
volume  and  surface  reverberation,  an  important  source  of 
interference  for  SAS  in  shallow  water  is  corruption  by  bottom 
reverberation  components  propagating  into  the  receiver  via 
surface-bounce  multipaths  (see  fig.  1).  Multipath 
interference  is  detrimental  to  SAS  image  quality.  Fig.  2 
shows  an  artificial  SAS  image  with  and  without  multipath 
interference  in  shallow  water  from  a  patch  of  reverberation 
containing  a  target  casting  a  shadow  on  the  seafloor.  The 
multipath  interference  adds  noise  to  the  SAS  image  and 
obscures  target  structures  such  as  shadows.  Another  effect 
is  the  inability  to  coherently  focus  the  multipath  components 
because  they  have  been  decorrelated  by  sea  surface 
roughness,  resulting  in  blurring  [2].  Multipath  interference 
becomes  particularly  severe  when  the  range  R  is  much 
greater  than  the  water  depth  L,  that  is,  R  »  L. 

In  principle,  a  vertical  array  receiver  could  be  used  to 
null  or  attenuate  the  multipath  interference  by  forming  delay- 
sum  beams  steered  toward  the  bottom  to  suppress 
interference  components  coming  from  surface  directions. 
Flowever,  the  closeness  of  the  multipath  angles  at  longer 
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Fig.  1.  Depiction  of  bottom  reverberation  components  propagating 
to  the  sonar  receiver  over  surface  bounce  multipaths. 
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Fig.  2.  Artificial  SAS  image  with  and  without  multipath  interference 
containing  a  target  that  casts  a  shadow  on  the  bottom.  The 
multipath  interference  fills  in  the  target  shadow. 

ranges  and  the  presence  of  higher-order  multipaths,  i.e., 
paths  involving  a  surface  and  bottom  bounce,  necessitates 
the  use  of  unrealistically  large  vertical  arrays  (space 
constraints  of  the  SAS  vehicle  limit  the  size  of  the  vertical 
aperture)  to  obtain  the  angular  resolution  needed  for 
sufficient  multipath  attenuation. 

Alternatively,  a  small  vertical  array  could  be  combined 
with  an  adaptive  method  such  as  a  generalized  sidelobe 
canceller  (GSC)  [3]  to  null  the  multipath  interference. 
Flowever,  application  of  adaptive  interference  cancellers 
here  is  complicated  by  the  simultaneous  presence  of  both 
signal  and  multipath  interference  (i.e.  signal  contaminated 
training  data)  which  results  in  slow  convergence  and  signal 
cancellation  (see  Boroson  [4]).  GSCs  are  also  prone  to 
signal  cancellation  from  errors  in  the  signal  blocking  matrix 
caused  by  array  calibration  and  signal  beam  pointing  errors. 
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Motivated  by  these  problems,  we  had  proposed  in  [5]  a 
method  for  blindly  separating  the  signal  from  multipath 
interference  that  estimated  the  separation  filter  weights  by 
maximizing  the  ping-to-ping  phase  coherence,  exploiting  the 
decorrelated  nature  of  the  multipath  interference,  rather  than 
a  least-squares  criterion  vis-a-vis  a  GSC.  It  works  well  while 
not  needing  precise  array  calibration  information  or  signal- 
free  training  data.  However,  the  preliminary  formulation  in  [5] 
required  multiple  pings  from  fixed  sonar  positions  for 
adaptation,  which  may  be  difficult  to  acquire  in  an  actual 
SAS  system  because  of  platform  motion.  In  this  paper  we 
extend  the  separation  methodology  to  a  moving  SAS 
system.  An  outline  of  the  paper  is  presented  below: 

2)  Problem  formulation:  The  mathematical 
representations  of  the  received  signal  plus 
multipath  interference  are  developed  and  we 
describe  how  pairs  of  vertical  array  measurements 
with  correlated  bottom  signals  can  be  obtained  by 
adjusting  the  forward  motion  of  the  SAS  vehicle  so 
that  the  array  phase  centers  overlap  between  any 
two  successive  pings. 

3)  Optimum  filtering  and  canonical  correlation- 
based  blind  separation:  We  review  optimum 
filtering  and  show  that  the  optimum  filter  weights 
can  also  be  determined  blindly  (without  needing  a 
signal  replica)  by  designing  filter  weights  to 
maximize  the  output  correlation  between 
overlapping  array  phase  centers,  taking  advantage 
of  the  correlated  nature  of  the  bottom  signal 
between  any  two  successive  pings. 

4)  Blind  estimation  of  filter  weights:  Instead  of 
using  canonical  correlation  analysis,  we  advocate 
estimating  the  filter  weights  by  minimizing  the 
circular  variance  of  the  phase  differences  between 
overlapping  array  phase  centers  and  motivate  why 
the  circular  variance  is  a  good  measure  of  phase 
coherence. 

5)  Numerical  Results:  Numerical  results  are 
presented  comparing  the  new  method  against  the 
CRLBs,  a  high  resolution  angle  of  arrival  estimation 
algorithm,  conventional  beamforming,  and  ideal 
null-steering,  showing  that  the  blind  method  works 
well  and  does  not  incur  much  performance  loss 
relative  to  parametric  methods. 

II.  PROBLEM  FORMULATION 

Suppose  the  SAS  system  has  a  planar  array  as  a 
receiver  and  a  point  source  transmitter.  For  the  sake  of 
simplicity,  (without  loss  in  generality)  assume  that  it  consists 
of  only  two  columns  of  equi-spaced  hydrophones  and  a 
transmitter  located  at  the  center  (see  fig.  3).  Denote  the 
complex  vertical  aperture  sensor  or  beam  outputs  (at  a 
given  frequency)  containing  the  backscattered  bottom  signal 


and  multipath  components  impinging  on  the  leftmost  and 
rightmost  staves  from  a  given  range  cell  in  the  kth  ping  cycle 
as  xk  and  yk  respectively.  We  can  then  model  the  stave 

p 

+  Z  cjJj + >  (2-1) 

7=1 

multipaths 
P  _ 

+  Z  c^j  +  nl ,  (2.2) 

J= i 

V 

multipaths 

where  the  vector  s  is  the  signal  basis  vector  or  array 
manifold,  hj  is  the  basis  vector  of  the  jth  multipath 
interference  component,  P  is  the  number  of  interfering 
multipaths,  ak  ,  ayk  are  the  complex  signal  amplitudes,  the 
cjk  are  complex  random  variables  corresponding  to  the 
interference  amplitudes,  and  the  elements  of  the 
background  noise  vector  nk  ,  n k  correspond  to  volume  and 
surface  reverberation  and  ambient  noise. 


outputs  as 


xk=aks 


signal 


and 


yk  =  gtf. 

signal 


•  —  Receive  hydrophones 

•  —  Transmitter 

•  —  Fictitious  monostatic 

receivers  based  on 
hyd  raphon  e-tran  smitte  r 
phase  centers 


Fig.  3.  Planar  array  geometry  and  creation  of  fictitious  monstatic 
elements  using  the  phase  centers  of  the  receiver-transmitter  pairs. 
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7  Overlapping  phase  centers 


Fig.  4.  Overlapping  vertical  phase  centers  obtained  by  adjusting  the 
ping  repetition  rate  and  the  SAS  vehicles  forward  speed. 
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Although  each  receiver-transmitter  pair  has  a  bistatic 
geometry,  they  can  be  approximated  by  an  equivalent 
fictitious  monostatic  receiver-transmitter  located  at  their  mid¬ 
point  (see  fig.  3)  [6,7].  This  approximation  holds  as  long 
A2  /  4r  «  A  where  A  is  the  distance  between  the  receiver 
and  transmitter,  r  is  the  range,  and  X  is  the  signal 
wavelength  [7].  That  is,  the  backscattered  signal  measured 
from  the  seafloor  by  the  bistatic  receiver-transmitter  pair  is 
approximately  the  same  as  that  would  be  measured  by  a 
single  receiver-transmitter  located  at  their  phase  center.  If 
we  adjust  the  speed  of  the  SAS  vehicle  and  the  ping 
repetition  rate  such  that  the  array  is  displaced  by  distance 
d/2  between  successive  pings,  the  array  phase  centers  will 
overlap,  as  depicted  in  fig.  4.  Therefore  the  stave  receiver- 
transmitter  pairs  corresponding  to  the  overlapping  fictitious 
monostatic  elements  will  measure  the  same  bottom  signal. 
This  is  the  principle  that  is  used  by  the  method  of  displaced 
phase  centers  for  SAS  motion  compensation  [6,7]. 

Thus  for  any  pair  of  overlapping  fictitious  staves,  we 
have  axk  «  ak+l .  Because  of  sea  surface  roughness,  we 
will  also  assume  that  the  multipaths  are  uncorrelated  from 
ping-to-ping  and  path-wise.  A  reasonable  model  for  the 
signal  and  multipath  interference  amplitudes  and  the 
background  noise  is  that  they  are  zero-mean  complex 
Gaussian  distributed  [3]  with  covariance  matrices 

Rx=Ry  =  AsssH  +  £  ^jhjhj  +  cr2I  (2.3) 
M_ _ _ / 

*7 

and 

Rxy  =  E[xky?+i]  =  PAsSSH  ,  (2.4) 


that  in  our  problem,  pairs  of  signal  and  noise  measurements 
(xk,  yk+l)  from  overlapping  phase  centers  are  observed  in 
which  the  signal  is  highly  correlated  between  the  vector 
pairs,  but  the  interference  is  uncorrelated.  It  turns  out  that 
these  differences  in  correlation  properties  between  signal 
and  interference  can  be  used  to  blindly  determine  the 
optimum  filter  weights  without  needing  knowledge  of  the 
signal  replica  [9].  Using  the  results  of  [9]  here,  the  filter 
weights  that  maximize  the  output  correlation 
E[wH xkyk+lw]  between  two  overlapping^phase  centers 
subject  to  the  constraint  that  Var[wxk]  =  1  are 
proportional  to  the  optimum  filter  weights  w  =  RJls  and 
can  be  calculated  from  Rx  or  Ry  and  Rxy  by  solving  the 
eigen-equation 


(3-1) 


for  A  and  w .  This  method  of  calculating  the  filter  weights 
is  based  on  the  technique  of  canonical  correlation  analysis 
from  multivariate  statistics.  The  statistician  Hotelling 
developed  canonical  correlation  analysis  in  the  1930s  for 
exploring  dependencies  or  relationships  between  random 
vectors  [10]  and  today  is  a  standard  statistical  analysis  tool. 
The  approach  is  useful  because  it  is  does  not  need  signal- 
free  interference  measurements  (i.e.  works  with  signal 
contaminated  data)  or  a  precise  signal  replica. 

However,  in  our  case  both  Rx  and  Rxy  are  unknown.  It  is 
easy  to  convince  oneself  that  plugging  estimates  of  Rx  and 
Rxy  into  (3.1)  and  solving  for  w  does  not  generally 
maximize  the  true  output  cross-correlation 


c(w)  =  ■ 


l2)I/2a>*yn+1 12) 


2x1/2 


(3.2) 


where  Xj  is  the  power  of  the  jth  path  and  p  is  the  correlation 
coefficient  between  ak  and  a£+l  . 

III.  OPTIMUM  FILTERING  AND  CANONICAL 
CORRELATION-BASED  BLIND  SEPARATION 

Our  objective  is  to  estimate  the  signal  expansion 
coefficients  ak  in  the  xk  (2.1)  in  the  presence  of  multipath 
interference.  The  minimum  least-squares  error  estimate  of 
the  signal  expansion  coefficient  ak  is  given  by 
ak  =  w”xk  where  wo  =  RJls  /(sH RJls)  [8].  These 
filter  weights  also  maximize  the  output  signal-to-noise  ratio 
[8].  This  of  course  assumes  that  we  know  the  interference 
covariance  matrix  /?7  and  the  signal  replica  s  .  Suppose  for 
a  moment  that  s  is  not  accurately  known  because  of  array 
calibration  errors  and  uncertainty  in  the  signal  arrival 
direction.  Is  there  a  way  to  determine  the  optimum  filter 
weights  without  needing  s  ,  that  is,  blindly  find  wo  .  Recall 


This  suggests  that  an  estimate  based  on  (3.1)  is  likely  to  be 
poor.  An  additional  concern  is  that  canonical  correlations 
needs  high  stability  in  both  magnitude  and  phase  of  the 
transmitted  and  received  signal.  Magnitude  stability  in 
practice  is  more  difficult  to  achieve  than  phase  stability 
because  of  transmitter  electronics  fluctuations.  Thus  we 
advocate  designing  the  filter  weights  by  directly  maximizing 
the  phase  coherence  between  the  (xk,yk+1)  pairs  rather 
than  indirectly  using  estimated  covariances.  In  the  next 
section  we  present  an  iterative  approach  for  estimating  the 
optimum  filter  weights  that  is  based  on  minimizing  the 
circular  variance  of  the  phase  differences  only. 

IV.  BLIND  ESTIMATION  OF  FILTER  WEIGHTS 

We  want  to  estimate  the  optimum  filter  weights  from  K 
independent  pairs  of  signal  and  interference  measurements 
(xx ,  y2  ), . . . ,  (xK ,  yK+l )  by  maximizing  the  output  phase 
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coherence  measured  between  the  pairs  using  some  metric. 
One  way  of  quantitatively  measuring  the  phase  coherence 
between  wHx^  and  is  to  calculate  the  phase 

difference  Zwxk-Zwyk+t  where  the  operator  Z  • 
means  the  argument  of  the  complex  number,  i.e.,  angle  or 
phase.  Following  this  idea,  a  reasonable  approach  then 
might  be  to  estimate  the  filter  weights  by  minimizing  the 
mean-square  phase  difference  between  the  overlapping 
phase  centers,  i.e.,  solving 

^  Z(zw»xk-Zw"yk+,y  ■  (4-D 

W  * 

However,  numerical  evaluation  of  phases  is  hard  because 
their  cyclical  nature  (modulo  2n )  may  require  phase 
unwrapping.  This  difficulty  can  be  overcome  by  using  the 
circular  variance  [11,12].  The  circular  variance  (CV)  is 
defined  for  given  a  set  of  angle  measurements  <h,...,<|>N  as 


cv[<t>}  =  l-V^+Z 

l  N 

x  =  —  Xcos^) 

^  n= 1 

j  =  Sin^) 

™  n= 1 


Based  on  the  above  discussions,  we  propose  estimating 
w  by  minimizing  the  circular  variance  of  the  output  phase 
differences, 


cv[w]  =  1  - 


J_y ^  ei(^wHxk-ZwHyk+l) 

Kfl 


(4.4) 


with  w  constrained  to  be  non-zero  to  avoid  getting  a  trivial 
solution.  Fig.  5  shows  a  comparison  of  the  proposed  blind 
method  to  the  conventional  GSC. 


Fig.  5.  Comparison  of  the  blind  method  against  the  conventional 
GSC  where  B  is  the  signal-blocking  matrix. 


and  lies  in  the  interval  [0,  1],  with  one  corresponding  to 
maximum  variance.  It  is  easy  to  evaluate  numerically  and 
commonly  used  in  analyzing  the  dispersion  of  angular  data 
such  as  wind  and  ocean  current  direction,  fracture 
orientations,  and  animal  movements.  The  CV  is  also  related 
to  the  maximum  likelihood  (ML)  estimate  of  the 
concentration  parameter  of  the  von  Mises  distribution  [11] 

pVM(</>)  =  [27rI0(k)r1eKCOS^)  ,  {)<(/>  <2tt  ,  (4.3) 

through  the  relation  1-_CV=I1(k)/I0(k)  where  k  is  the 
concentration  parameter,^  is  the  mean  angle,  and  lm(k)  is 
the  modified  Bessel  function  of  order  m.  The  von  Mises 
distribution  is  often  used  to  model  circular  or  angular  data 
and  is  considered  to  be  the  circular  analog  of  the  regular 
Gaussian  distribution  [11,12]  and  approaches  the  Gaussian 
distribution  as  ft*  — »  oo  [11].  This  suggests  that  the  CV  is 
likely  a  good  measure  of  phase  dispersion  or  coherence. 
Viterbi  also  showed  that  (4.3)  is  the  probability  density  for 
the  carrier  phase  estimate  of  a  sinusoidal  signal  in  noise 
obtained  by  first-order  phase-locked  loop  and  that  the 
concentration  parameter  k  is  proportional  to  signal-to-noise 
ratio  [13].  The  CV  has  also  been  previously  used  in 
frequency  estimation  and  can  be  interpreted  as  a  type  of 
circular  regression  in  that  problem  at  high  signal-to-noise 
ratios  [14]. 


Since  no  closed  form  solution  is  available  for  the  above 
minimization  problem,  some  form  of  iterative  search  is 
required,  e.g.  steepest  descent,  quasi-Newton  etc  [15].  We 
chose  to  use  the  Fletcher-Reeves  (FR)  method  [15]  for 
performing  the  optimization.  The  search  direction  for  the  kth 
FR  iteratation  is  given  by 

dk=-gk+{\gkfFl\gk-xtF^k-x  -  (4-5) 

where  gj  is  the  gradient  of  the  phase  variance  cv[w] 
evaluated  using  the  iterate  from  the  jth  step.  The  updated 
weight  vector  estimate  wk  =  wk_x  +  Adk  is  found  by 
minimizing  cv[wk_x  +  Adk  ]  with  respect  to  A.  This 
method  is  simple  and  has  quadratic-like  convergence  [15]. 
We  will  discuss  starting  points  and  convergence  in  more 
detail  later. 

This  procedure  is  “blind”  because  the  signal  replica  s  is 
not  needed  for  estimating  the  weight  vector.  But  wean  be 
only  be  uniquely  determined  up  to  a  scale  factor.  If  the  time 
series  is  to  be  reconstructed  after  doing  the  separation  in  the 
frequency  domain,  then  at  least  some  partial  knowledge  of 
the  array  geometry  and  calibration  is  needed  to  maintain  unit 
gain  signal  response  for  the  signal  reconstruction. 
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A.  Optimality  of  the  Circular  Variance 

For  a  given  set  of  filter  weights,  formula  (4.4)  is  a 
function  of  the  phase  differences  between  two  correlated 
complex  Gaussian  random  variables,  i.e., 
(f)  =  ZwHxk  —ZwHyk+l.  Their  probability  density  is  given 
by  (see  [16]  and  also  [17]) 


pCG(</>\w) 


27T(\-P2f12 


where / 3  =  pcos((f)  —  (f))  and  p  is  the  correlation  between 
wHxk  and  wHyk+1.  By  inspection  (4.6)  is  not  the  same 
as  the  von  Mises  density  (4.3).  Since  we  are  implicitly 
seeking  to  maximize  the  correlation  between  the  filter 
outputs  by  minimizing  the  phase  circular  variance,  a 
question  arises  as  to  how  well  the  CV  characterizes  the 
correlation  between  complex  Gaussian  random  variates. 
Since  density  (4.6)  is  a  function  of  the  correlation  coefficient 
p,  it  would  be  desirable  prove  sufficiency  [18]  or  at  least 
approximate  sufficiency  of  the  CV  for  p.  However,  this 
appears  difficult.  In  appendix  A,  we  looked  at  the  related 
hypotheses  testing  problem  of  determining  whether  a  given 
population  is  distributed  as  p=0  (uniform)  or  p>0  using  1-CV 
and  the  likelihood  ratio  (LR)  as  test  statistics  and  show 
numerically  that  the  CV  appears  to  “capture”  much  of  the 
information  regarding  p. 


V.  NUMERICAL  RESULTS 


Here  we  present  some  numerical  results  investigating 
the  error  surface  and  comparing  the  blind  separation  method 
by  computer  simulation  against  high  resolution  angle 
estimation,  Cramer-Rao  lower  bounds  [18],  conventional 
beamforming  and  the  clairvoyant  null-steerer  and  an  artificial 
SAS  image  example.  We  begin  by  discussing  the 
experimental  scenario. 

The  SAS  system  is  modeled  as  having  a  10  element 
equi-spaced  vertical  array  receiver  and  ensonifies  the 
bottom  with  a  monochromatic  signal  pulse  in  each  ping 
cycle.  One  signal  is  directly  backscattered  from  the  bottom 
to  the  receiver  and  another  arrives  to  the  receiver  via  a 
surface  bounce  path,  with  both  arrivals  modeled  as  plane 
waves  (see  fig.  1).  Two  overlapping  phase  center  pairs 
(xk9yk+1)  are  collected  in  each  SAS  cycle  and  are 
modeled  as  being  jointly  complex  Gaussian  distributed  with 
the  covariance  and  cross-covariance  matrices  of  (2.3)  and 
(2.4)  respectively.  Their  parameters  were  set  to  X.s=1 ,  p=.95 
(ping-to-ping  signal  correlation),  7^=1  (the  signal  and 
multipath  interference  have  the  same  power),  and  cr2=.01. 
The  new  blind  method  was  implemented  in  beamspace 
analogous  to  that  of  a  GSC,  i.e.,  using  the  outputs  of  two 
beams,  one  steered  at  the  desired  bottom  signal,  and  a 
interference  reference  beam  orthogonal  to  the  signal  steered 


at  the  multipath  interference  (like  in  fig.  5).  This  is  equivalent 
to  operating  on  the  transformed  data  uk  —BHxk  and 
vk  =  BH  yk+1 ,  where  B  is  a  10  x  2  matrix  whose  columns 
correspond  to  the  bottom  and  surface  steering  replica 
vectors. 

A.  Circular  Variance  Surface 

In  the  first  experiment,  we  examined  the  shape  of  the 
circular  variance  surface  as  a  function  of  the  filter  weights  for 
insight  into  the  difficulty  of  using  gradient-based  iterative 
search  methods.  Since  only  one  multipath  component  is 
present,  the  weight  vector  in  formula  (4.4)  can  be  written  as 
w  =  [1  (ar  +  iat)]T  .  K=10,  20,  30,  and  40  independent 
data  snapshots  pairs  were  computer  generated  using  (2.1) 
and  (2.2)  and  plugged  into  the  circular  variance  formula  (4.4) 
and  plotted  in  fig.  6  as  a  function  of  aR  and  ah  We  see  that 
as  the  number  of  data 


K=1Q  K=20 


Fig.  6.  The  output  circular  variance  plotted  as  a  function  of  aR  and  sl\ 
for  K=10,  20,  30  and  40  independent  data  snapshot  pairs.  The 
symbol  “*”  corresponds  to  the  optimum  filter  weight  values. 

snapshots  increases,  the  surface  tends  to  become  more 
regular  or  quadratic-like  with  the  circular  variance  minima 
being  close  to  that  of  the  optimum  filter  weights.  However, 
surface  irregularities  and  local  stationary  points  do 
occasionally  occur  outside  the  concave  region  surrounding 
the  global  minima  as  seen  in  fig.  6,  suggesting  that  it  is 
important  to  have  a  good  starting  weight  vector  to  avoid 
getting  stuck  in  local  stationary  points  during  the  iterative 
search.  Some  possible  approaches  for  obtaining  good 
starting  vectors  for  the  Fleetcher-Reeves  method  [15] 
include  using  the  GSC  filter  weight  estimate  or  its  gradient 
as  an  initial  search  direction  or  using  prior  information  about 
the  geometry  (sonar  and  water  depth,  and  range)  or 
computing  the  differences  xk  -  yk+l  to  obtain  a  nearly 
signal-free  data  set  to  estimate  an  initial  guess  (it  can  be 
argued  that  if  the  signal  is  highly  correlated  between  the 
phase  center  pairs,  this  procedure  will  remove  most  of  the 
signal  assuming  no  transmitter  magnitude  fluctuations). 
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B.  1 -Multipath  Performance  Evaluation 

The  performance  of  the  blind  method  is  now  compared 
against  the  Cramer-Rao  lower  bound  (CRLB)  (see  the  CRLB 
derivation  in  appendix  B),  high  resolution  angle  of  arrival 
estimation  algorithm  of  Kumaresan  and  Tufts  [19],  the 
clairvoyant  null-steerer  (constructed  using  perfect 
knowledge  of  signal  and  interference  arrival  directions),  and 
the  conventional  beamformer.  30  independent  data  pairs 
were  computer  generated  in  each  trial  according  to  the 
experimental  scenario  described  at  the  beginning  of  this 
section  and  both  the  blind  method  and  the  minimum-norm 
methods  were  implemented  (the  minimum-norm  method 
estimated  the  signal  by  first  estimating  the  signal  and 
interference  spatial  frequencies  and  then  fitting  plane-waves 
of  those  frequencies  in  a  least-squares  sense  to  the 
observed  data).  The  blind  method  was  implemented  using 
an  initial  guess  of  w0  =  [1  0]T  to  start  the  Fletcher- 
Reeves  algorithm.  Both  the  error  in  estimating  the  weight 
vector  (nominal  weight  vector  subtracted  frqm  estimated 
weight  vector)  and  the  signal  (a*  =wHxl)  were 
tabulated.  A  total  of  two  thousand  independent  computer 
simulation  trials  were  performed  at  different  spatial 
separations  to  measure  the  mean-square-errors  (MSEs). 


Fig.  7.  The  measured  MSE  of  the  blind  method  (weight  vector  and 
signal  estimation  errors)  plotted  as  a  function  of  signal  and 
interference  separation  in  DFT  bins  and  compared  against  the  high 
resolution  angle  of  arrival  estimation  algorithm,  the  CRLB, 
conventional  beamformer,  and  clairvoyant  null-steerer. 

The  lowermost  plot  of  fig.  7  compares  the  signal 
estimation  MSE  of  the  blind  method  against  the  high 
resolution  method,  the  CRLB  (see  appendix),  clairvoyant 
null-steerer,  and  conventional  beamformer  as  a  function  of 
signal  and  interference  vertical  separation  in  DFT  bins.  We 
see  that  the  blind  method  performs  well  with  respect  to  the 
other  methods,  and  becomes  very  competitive  against  the 


parametric  high  resolution  method  once  the  signal  and 
interference  are  within  a  1/2  DFT  bin  width  of  each  other.  We 
also  see  that  the  blind  method  signal  MSE  is  fairly  close  to 
that  predicted  by  the  CRLB  (see  appendix  B).  The 
uppermost  plot  of  Fig.  7  compares  the  weight  vector 
estimation  error  against  the  CRLB  for  the  weight  vector  (see 
appendix  B).  Although  the  blind  method  does  not  approach 
the  CRLB,  we  must  bear  in  mind  that  the  CRLBs  were 
derived  assuming  that  both  signal  magnitude  and  phase 
information  are  used  in  the  estimation  process,  whereas  the 
blind  method  uses  only  phase  information.  Next,  the 
experiment  was  repeated  with  the  spatial  separation  set  to 
1/2  DFT  bin  and  the  number  of  data  snapshots  pairs  used  for 
estimating  the  weight  vector  varied  from  10  to  100  (see  fig. 
8).  As  expected,  the  uppermost  plot  of  Fig.  8  shows  that 
blind  method  does  not  approach  the  weight  vector  CRLB  as 
the  number  of  snaphot  pairs  K  increases.  Flowever,  in  terms 
of  the  signal  estimation  MSE,  the  blind  approach  performs 
well  and  is  close  to  that  of  the  high  resolution  angle 
estimation  method. 


Fig.  8.  The  measured  MSE  of  the  blind  method  (weight  vector  and 
signal  estimation  errors)  plotted  as  a  function  of  the  number  of  data 
snapshot  pairs  for  used  adaptation  and  compared  against  the 
CRLBs  and  the  other  methods. 

C.  2-Multipath  Artificial  SAS  Image  Example 

In  the  last  example,  we  simulate  a  two  multipath 
interference  situation  (a  surface  and  surface-bottom  path) 
and  show  that  the  blind  method  outperforms  the  generalized 
sidelobe  canceller  and  the  conventional  beamformer  in 
terms  of  SAS  image  quality.  The  same  experimental 
scenario  is  used  as  before,  except  for  the  presence  of  a 
surface-bottom  path  and  an  object  casting  a  shadow  on  the 
seafloor.  The  direct  path  bottom  signal  was  separated 
vertically  from  the  surface  path  in  spatial  frequency  by  .4 
DFT  bins  and  by  -.7  DFT  bins  from  the  surface-bottom  path. 
The  surface  and  surface-bottom  path  reverberation  levels 
were  set  to  the  same  power  and  a  total  of  40  data  snapshot 
pairs  were  collected  for  processing  while  the  SAS  system 
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moved  on  a  linear  track.  The  GSC  was  implemented 
erroneously  (signal  reference  and  signal-blocking 
interference  reference  beams  constructed  using  incorrect 
signal  direction  information)  with  the  assumed  signal 
direction  mismatched  by  .2  DFT  bins  away  from  the  actual 
signal  direction. 

Fig.  9  shows  the  SAS  images  generated  using  the 
conventional  beamformer  (a  beam  steered  toward  the 
bottom),  the  blind  separation  method,  and  the  GSC.  Note 
that  the  object  shadow  occurs  in  range  cells  30-40.  By  visual 
inspection,  we  see  that  the  object  shadow  in  the 
conventional  beamformer  image  is  barely  perceptible,  i.e., 
obscured  by  the  multipath  interference.  The  GSC  image 
shadow  quality  is  somewhat  better,  although  the  shadow 
area  still  contains  considerable  multipath  leakage.  The  SAS 
image  obtained  from  the  blind  method  has  the  best  shadow 
quality  out  of  the  three  methods.  The  GSC  performed  poorly 
because  of  signal  contamination  in  both  the  signal  reference 
beam  and  the  signal-blocking  interference  reference  beams. 


range  cell  number 


APPENDIX  A:  OPTIMALITY  OF  THE  CV 

One  approach  for  getting  insight  into  the  sufficiency  of  the 
CV  for  p  is  to  look  at  the  hypotheses  testing  problem  of 
determining  whether  K  independent  phase  difference 
measurements  are  distributed  according  to  p=0  (uniform)  or 
p>0  using  1-CV  and  the  likelihood  ratio  (LR)  (constructed 
using  density  (4.6)  and  assuming  the  mean  phase  is 
unknown)  as  test  statistics.  We  implemented  and  evaluated 
the  performance  of  the  two  detectors  for  5,  10,  20,  and  30 
independent  and  identically  distributed  phase 
measurements  (number  of  measurements  is  denoted  as  K  in 
the  paper)  using  5000  independent  computer  simulation 
trials  and  tabulated  the  probability  of  detection  as  a  function 
of  correlation  p  with  the  probability  of  false  alarm  set  to  10'2. 
The  results  in  fig.  10  show  that  the  CV  performs  nearly  as 
well  as  the  optimum  LR  detector,  suggesting  that  the  CV  is  a 
good  measure  of  correlation  for  this  problem  and  is  likely 
nearly  sufficient  at  least  in  these  examples. 


Fig.  9.  Artificial  SAS  images  created  after  multipath  rejection  using 
conventional  beamforming  (top),  blind  separation  (middle),  and 
GSC  (bottom).  Note  that  all  of  the  methods  were  applied  in  the 
vertical  domain  and  the  true  object  shadow  occurs  between  range 
cells  30-40. 


Fig.  10.  The  measured  probability  of  detection  vs.  p.  The  blue  and 
red  dashed  lines  correspond  to  the  likelihood  ratio  and  circular 
variance  test  statistics  respectively. 


APPENDIX  B:  CRAMER-RAO  LOWER  BOUNDS 


VI.  CONCLUDING  REMARKS 

These  results  suggest  that  blind  separation  of  the  signal 
from  the  multipath  interference  is  competitive  with  high 
resolution  angle  of  arrival  estimation  algorithms  and 
compares  favorably  to  the  CRLB-predicted  signal  MSE.  An 
important  advantage  of  doing  the  separation  blindly  is  that 
signal-free  training  data  and  highly  accurate  array  calibration 
information  are  not  needed,  unlike  high  resolution  methods 
or  GSCs.  Flowever,  more  work  is  needed  in  understanding 
the  convergence  properties  of  the  iterative  scheme  and 
performance  as  the  number  of  multipaths  increases. 


Flere  we  describe  the  steps  for  deriving  the  Cramer-Rao 
lower  bounds  on  the  error  in  estimating  the  filter  weights  and 
bounding  the  signal  estimation  error  from  K  independent 
pairs  of  signal  and  interference  measurements 
(x1? (x^,j^+1)  for  the  two  channel  case.  We 
assume  that  nothing  is  known  about  the  signal  and 
interference.  Suppose  for  a  moment  that  there  is  no  additive 
white  noise  (a2=0)  and  the  signal  is  perfectly  correlated 
between  the  pairs  ( xk ,  yk+l) .  Then  the  filter  weights  which 
provide  the  best  phase  coherence  (or  set  the  CV  to  zero) 
must  be  the  null-steerer 
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The  optimum  filter  weights  wo  =  RJls  also  approach  (B.1) 
as  the  multipath  interference  becomes  much  stronger  than 
the  background  noise,  i.e.,  \  »  a1  [20].  From  (B.1), 
another  interpretation  of  the  weight  vector  estimation 
problem  is  one  of  estimating  the  part  of  the  signal  h±  that 
lies  in  the  orthogonal  complement  of  the  space  spanned  by 
the  interference  basis  vectors. 

_ln  the  two  channel  case,  if  h=[-a  l]r,  then  we  have 
h±=[la]T.  Using  the  above  interference-based 
coordinate  system,  the  signal  basis  vector  can  now  be 
expanded  as 

s  =  (hL  +  J3h)/ 

Recall  that  the  joint  complex  Gaussian  distribution  of  the 
data  pairs  (xk,yk+ j)  is  characterized  by  the  covariance 
matrices  of  (2.3)  and  (2.4).  Plugging  (B.2)  into  (2.3)  and 
(2.4)  (a2  cannot  be  estimated  uniquely  here  because  we 
only  have  two  channels  or  degrees  of  freedom),  the 
unknown  parameters  to  be  estimated  are  a,  [3,  Xs,  K  and  p. 
The  elements  of  the  Fisher  information  are  easy  to  evaluate 
for  this  problem  and  are  given  by  [18] 


/i_i_  +  j3h 
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where  R  =  cov([xf  I  yl+l]T )  is  the  nominal  data 
covariance  matrix  and  the  notation  0m  means  mth 
parameter.  The  bound  on  the  mean  square-error  in 
estimating  the  mth  parameter  then  is 

y‘L  (B.4) 


It  is  straightforward  to  evaluate  (B.3)  and  (B.4)  numerically 
to  obtain  the  bound  on  the  error  in  estimating  the  weight 
vector  w  or  its  parameter  a. 

We  can  also  use  the  weight  vector  CRLB  to  glean  insight 
into  the  signal  estimation  error.  If  the  interference  covariance 
matrix  is  R|  and  the  bound  in  estimating  a  is  then 
approximately 

£[Ih^  -ak  I2]  »(l+l  <%  \2)(f  +o2\piRIB\12  (B.5) 


where  a0  is  nominal  value. 
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